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Abstract 

We discuss the excitation of sterile neutrinos in the early universe using general 
quantum kinetic equations, which also incorporate a possible primordial mag- 
netic field B. We find a new contribution to the excitation propability, which 
has its origin in the shrinkage of the spin vector. In the absence of B nucle- 
osynthesis implies the constraint | Am 2 | sin 2 26q < 1.6 x 1CP 6 eV 2 which is 
more restrictive than previous estimates. We also present examples of possible 
stringent limits for B ^ 0. 
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1. Introduction. 



The recent observation of the cosmic microwave background temperature anisotropy 
on large scales by the COBE satellite has hinted towards the existence of a hot 
component in the dark matter (HDM) of the universe A light neutrino is an 

obvious candidate for hot dark matter. As is well known, there are a number of unre- 
lated astrophysical and cosmological constraints on light neutrinos, which have to do 
with solar neutrinos and the deficit of muon neutrinos in the atmospheric neutrino 
fluxes ||. Reconciling all these issues in the context of a three generation mixing 
model only is quite difficult, and as was first pointed out in ||, if neutrinos are to 
solve the dark matter, solar and atmospheric problems simultaneously, one actually 
has to introduce a fourth neutrino, which could well be a sterile neutrino. Moreover, 
it has been suggested || that the recently observed signals in the KARMEN detector 
could be explained by the the decay pattern 7r + — > /i + zA,, where v s is light sterile 
neutrino. Simple extensions of the standard electroweak model which can accommo- 
date all known hints for neutrino masses, including solar and atmospheric neutrino 
observations also postulate the existence of a light neutrino v s |5|||. In some of these 
models such sterile neutrino is the HDM candidate M. (Sterile neutrino as warm 
dark matter has recently been discussed in 

Very tight constraints on the neutrino mass matrix that includes a singlet v s can 
be deduced from the primordial nucleosynthesis bound on the excess relic energy 
density at the proton- neutron freeze-out, which takes place at T w 0.7 MeV. The 
bound derives from the evaluations of the abundances of primordial 4 He and D and 
is usually quantified in units of relativistic neutrino species. For a recent discussion 
on the neutrino bound, see |§ (see also |§). 

Because of the mixing, sterile neutrinos can be produced in weak collisions and 
thus be brought into equilibrium, whence they would count as an extra relativistic 
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neutrino degree of freedom, in clear contradiction with the observations | 10| . Equi- 
libration can be avoided for a certain range of the mixing parameters, though. The 
region of the active-sterile neutrino oscillation parameters Am 2 = m 2 — m 2 and sin 2 29 



excluded by the Big Bang Nucleosynthesis (BBN) has been estimated to be |TT 

| Am 2 | sin 4 29 > 5 x KT 6 eV 2 , v a = u e , 
| Am 2 | sin 4 29 > 3 x 1(T 6 eV 2 , u a = v^ T . 



Folding in an actual nucleosynthesis code into neutrino evolution equations has ver- 
ified this result |12| . 



Strong random magnetic fields change neutrino oscillations drastically. This is a 
pertinent issue because of the possibility that the observed galactic magnetic fields, 



which are of the order of microgauss, might have a primordial origin JT3| . If this 
is true, then the plasma of the early universe sustained enormous magnetic fields 
which may have affected also particle interactions. It is natural to assume that the 
primordial field is random. This is because any magnetic field is imprinted on the 
almost chaotically comoving plasma, which entangles and mixes the field lines. Re- 
cently it was pointed out fl4| that in presence of such random magnetic fields the 
BBN constraints on sterile neutrino become more stringent than in isotropic plasma. 

In the present paper we investigate the kinetics of the active-sterile neutrino con- 
versions in hot plasma of the early universe using a set of quantum kinetic equations. 
(For earlier applications, see |p!0HlTip!5| .) Our formalism is general and takes also into 
account quantum damping as well as the effects of randomly magnetized plasma. 
Apart from reproducing the earlier results obtained in particle approach, the kinetic 
approach reveals new effects, which give rise to a novel constraint on the neutrino 
mixing parameters. It turns out to be more stringent than the previous ones, which 
have been obtain ed from one-particle Schrodinger equation. 

The organization of the paper is as follows. To describe neutrino propagation 
in hot plasma we derive, in Section 2, neutrino dispersion relations in magnetized 



plasma. We also discuss briefly the behaviour of random magnetic fields. In Section 
3, starting from the quantum kinetic equations (QKE) first derived in ||16|| , we obtain 



a general analytic solution for the neutrino conversion probability in equilibrium 
plasma at temperatures T 3> m e . Special care is taken to include correctly the 
effects of a possible random backround magnetic field. In Section 4 we apply our 
general formalism to the special case of isotropic plasma. We find a new shrinking, 
aperiodic analytic solution for u e — > u s neutrino conversions, which gives rise to a new, 
stringent constraint on the mixing parameters. In Section 5 we include a primordial 
magnetic field and find for large fields a simple analytic expression for the conversion 
probability. The solution is a generalization of the one-particle result of [fPjl . 

In Section 6 we discuss the results and their applications both to neutrino physics 
and to the astrophysics of galactic magnetic fields. 

2. Neutrino propagation in hot plasma 
2.1. Random magnetic fields 

In order to describe active-sterile neutrino oscillations in a general medium, which 
may also contain a random magnetic field, we need neutrino dispersion relations. We 
also need a model for the primordial magnetic field. Although there are a number of 
suggestions as to how a large magnetic field could arise in the early universe ||l7| , |T8f , 
we prefer a more phenomenological approach. 

A primordial magnetic field, no matter what its origin, is imprinted on the comov- 
ing plasma, which in the early universe has a very large conductivity. We shall assume 
that the primordial plasma consist of elementary magnetic domains of size L . Within 
each such domain the magnetic field is taken to be uniform and constant, and the 
field in different domains is randomly aligned. For the root mean squared magnetic 
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field B rms = sj (B 2 ), averaged over a volume L 3 3> Lq, we assume the scaling law 

ft™ = *(£)'(£)'. (2 ,, 

The temperature dependence of B rms reflects simply magnetic flux conservation. 
How does B rms scale with distance is an unsolved issue and reflects our ignorance 
e.g. on the question as to how uncorrelated the fluxes in the neighbouring magnetic 
domains actually are. In the case of no correlation statistical averaging gives for the 
parameter p the value p = 1/2, whereas a random- walk argument yields the scaling 
P = 3/2 H. 



Dissipation of B was investigated in [1£,2C]. There it was found that the dis- 
sipation length at the recombination time is about 10 10 cm, which corresponds to 
10 10 cm x (Tj.ee/Te_Bv) — 2 x 10 4 cm at nucleosynthesis. If primordial magnetic 
fields are to be the seed field for the galactic dynamo [0 , the primordial field should 
survive until recombination. This sets a limit for the size of the random magnetic 
field domain L : 

L > L™ m = 10 3 cm x (^) . (2.2) 



Primordial nucleosynthesis may be used to set limits on the size of B fl20|,|2T 



For an uncorrelated random magnetic field the mean field (Bj) = 0, whereas 

(5 l (x)5 J (x')) = ^M (3) (x-xO. (2.3) 

Here the correlation length A is determined by the domain size L and the value of 
root mean squared field B rms at the horizon scale L = Ih(T) |p~9|] : 

1 3 

t = —[7, — 7rrB rms (l H )L , p ^ 3/2 , 
A vr(3 - 2p) 



\ = -{\n l -f)B 2 rms {lu)Ll, p = 3/2 
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The horizon scale defines a cut-off k max for the wave number. Physically it corresponds 
to the scale of the inhomogeneity, k max = 27r/L ||19|| . This is the reason for the 
appearance of the horizon size in the above formulas. 

2.2. Neutrino dispersion relations 

The ultrarelativistic dispersion of a standard electroweak neutrino in hot plasma 
depends on the neutrino interaction potential, which consists of two parts: 

V = V^ vec ^ + y( axial \ (2.5) 

The potential is determined by the neutrino forward scattering amplitude off all 
particle species in the plasma, including magnetized charged leptons and antileptons. 
The vector interaction potential V^ vec ^ for a neutrino with a momentum p = (p) = 
3.15 T is given by [p2j 



V wc = G F V2r h [n aaym -A 1 



where n asym depends on the particle-antiparticle asymmetries in the plasma, normal- 
ized to the photon density n 7 = 0.244T 3 , and A ~ 55. In the hot primordial plasma 
with T ^> m e , where particle asymmetries are small, the vector interaction potential 
Kec is dominated by the second, non-local term, which reads 

IV^I ~ 3.4 x l(r 20 (^^) 5 MeV . (2.6) 

The axial potential Kxial is present only if the plasma supports a magnetic field. It 
is given by [14| 

v (axicd) _ „ k - B , Vlff ( R 2 Q- B ) 2 \ (0 7 n 
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Here the quantity j2 e ff is defined by 

n eff = — ^ AJ ~ -12 c A x 10- 13 n B (—— 

p " v^vr 2 VMeV 

where ca = TO. 5 is the axial coupling in the weak lepton current (the upper sign is 
for v e , the lower one for v^ T )-> and /is = e/2m e is the Bohr magneton. In our case 
it is sufficient to consider only the first term of y( axial ) since according to eq. (ZTT) 
HeffB <C ~ 3T. 

3. Neutrino conversions in hot plasma 
3.1. Averaging over random magnetic fields 

Quantum kinetic equations take into account the inherent quantum nature of neu- 
trino oscillations ( [16| and references therein, see also ||2~3||). Adopting this approach 



we now derive a general expression for the probability of active-sterile neutrino con- 
versions in the early universe. 

The time evolution of the system of the active neutrino v e and a sterile neutrino v s 
can be described in terms of a polarization vector (P (t), P(t)), whose z-component 
gives the the excess of v s over v e in the neutrino ensemble at a given moment of 
time. In Ref. [^6] the general time evolution equations of the polarization vector 
were derived (eqs. (28) to (31) of [|T6l] ) , and these are applicable for any neutrino 
transitions. In our case, where one of the neutrinos is sterile with no interaction with 
the background matter, these equations are considerably simplified. 

We shall assume a hierarchy between the relaxation times of the spatial distri- 
bution along the neutrino trajectory, and the momentum distributions. This as- 
sumption allows us to factorize out the equilibrium momentum distribution f u (k) = 
[exp(fc/T) + l] -1 , so that the density matrix can be written as (k is the unit vector 



6 



in the direction of k) 



p„(f,k) = \[P (t) + PifakWMk) . (3.1) 



The distributions of the active neutrino u e and the sterile neutrino v s can be written as 

f Ue (t,k) = hp (t) + p z (t 1 k)]f u (k) , 

(3.2) 

f Us (t,k)= 1 -[P (t)-P z (t,k)]fu(k) • 
Let us note that our definitions ( |3.1|) and ( p.2|) are more general than those given 
in 



16] 



There is a hierarchy of scales in our problem: Lq <C Lw "C Ih, where Lq is the 
previously defined scale of elementary domains in the randomly distributed magnetic 
field, Lw is the length scale of weak interactions and Ih is the horizon scale. There- 
fore we should average our QKE over the random magnetic field distribution before 
integration over the momenta. Moreover, we are interested in the neutrino conver- 
sion probability at physical distance much larger than L . It is described by the 
mean value of the flavour space polarization vector ^P z (t, k)^ = P z (t), which does 
not depend on the neutrino momentum. 

As the temperature changes slowly as compared with the neutrino oscillation rate 
one can write, using factorization (|3.2j), the active-sterile evolution equations in the 
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dP (t) 
dt 

dP x (t,k) 
dt 

dP y {t,k) 
dt 

dP z (t,k) 



R(t,k)/f v (k) , (3.3) 

-V z (t,k)P y (t,k) - D(t,k)P x (t,k) , (3.4) 

V z (t, k)P x (t, k) - P z (t, k)V x - D(t, k)P y (t, k) , (3.5) 

df V x Py(t,k)+R(k,t)/f v (k) . (3.6) 

Here R(k,t) is the rate of the annihilation processes {yj? e <-» e + e~), the general 



form of which is given in [lq|. The damping parameter D(k,t) is determined by the 



inelastic- and elastic neutrino collisions. 



The neutrino interaction potential V z is given by 

V e (t, k) = V^' ec) - A cos 26 + Mf B \\ (*) ( 3 - 7 ) 

where V vec is given by the eq. ( |2.6| ), the magnetic term is approximated from eq. 
(13), A = (m\ - ml) /2k, and V x = Asin2# . 

The system of eqs. (|3.3|) yields the following quite complicated integro-differential 
equation for P z : 



[P z (t) + V x 2 P z (t) + (P z (t)- 



R(t). 



D(t)]F(t) 



+ {y^ - A cos 20, 



(! ) / dti 

o 



vdP z {t x ) R(h 



dti f, 
dP z (h) R(h) 



F{h) 



dt\ 



fu 



+<yW - Aoos26lo)AW/{B,|(t) ^^(^^ - ^-)F{ti) 







(3.8) 



where Fit) = exp(/ * D(t' )dt') (here we omit in all the arguments the momentum 
variable k.) In the case of a random magnetic field integration over the dimen- 
sions transversal to the neutrino pragation direction allows us to omit the terms 



linear in B». For the product -B||(t).B||(ti) we may use the approximate result [14 
B\\(t)B\\(t\)^ I (.B^\ ~ L 5(t — ti), so that this term reduces to 

dP z (t) 



2Y 



R(t)/f v exp( / D(t 2 )dt 2 ) . 



dt J Jo 

Here T is the damping induced by the presence of a magnetic field. Using eq. ( [2T4] 
for p 3/2 yields the following expression: 

3 



r = K// ( B \\ L ° 



(3.9) 



4tt(3 - 2p) 

Differentiating eq. (|3.8| ) with respect to time, and omitting the common exponen 
tial factor, we finally obtain the following third order differential equation: 



S 



/(*) 



d 3 p z 



dt 3 



+ 2(r + D(t, k)) 



d 2 P z 
dt 2 



+ [A 2 m + 2TD(t,k) + D 2 (t, k) + 



+2T 



R(t,k)[2TD(t,k) + D 2 (t,k) + 
dR(t, k) 



dD(t,k)\ dP z 
dt J dt 
dD(t,k) 



D(t,k)V 2 P z 



dt 



+ (V {vec) - Acos29 f 



dt 



where A m is the usual neutrino oscillation frequency in isotropic medium: 



A m = [(0^-Acos2^o) 2 + K 2 ] 1 



/2 



(3.10) 



(3.11) 



3.2. Master equation for active-sterile conversion 

Before neutrino decoupling, pair annihilation does not contribute to the factor R, 
but it contributes, like elastic neutrino collisions, to the collision damping coefficient 
D(t, k). At fixed temperature D does not depend on time. This allows us to rewrite 
the evolution equation of the active-sterile neutrino conversion probability 



P = P» 



in the following simplified form: 



1 - 



Po(0) 



^ + 2(r + rv)^ + (Am + 2rrv + r 2 m )^ + t w v 2 p 



t w v 2 



(3.12) 



(3.13) 



Here we have averaged eq. ( |3.10| ) over the normalized momentum distribution, and 
used the notations 



w 



V 2 
1 wi 



dkf(k)D(k) = (D), 
J dkf(k)D 2 (k) = (D 2 ) 



(3.14) 



where dk = nd 3 k/(2n) 3 . 
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Eq. ( |3.13|) is our master equation for the active-sterile neutrino conversion. The 



initial conditions are easily determined from the equations (|3.3|)-(|3~3]) and are given by 



P(0)=0, P(0)=0, P(0) = ^f. (3.15) 



The solution of eq. (|3.13|) is of the general form 



P(t) = - + C ie klt + C 2 e k2t + C 3 e M . (3.16) 



Here ^1,2,3 are the roots of the Cartan equation 

k 3 + a 2 k 2 + a%k + clq = 
where the constant coefficients read 



a 2 = 2(r + rv) , 

a 1 = A^ l + 2rr w + r^ 1 , (3.17) 
a-o = ^wVx ■ 



The roots 



2{Y + T W ) 
h = - h (si + s 2 ) , 

k 2 = ~ ^{si + s 2 ) + —(s 1 -s 2 }, (3.18) 

h = _^£±I^1 _ I( Sl + S2 ) _ % -^-{si - s 2 ) , 

can be found by solving the Cartan equation in the standard way so that s\ = 
(r + ^fr 2 + g 3 ) 1 / 3 , s 2 = (r — y/r 2 + g 3 ) 1 / 3 with r = (aia 2 — 3a )/6 — a 2 /27, q = 
ai/3 — a|/9. The constants Cj = Ci(ki, k 2 , /c 3 ) in eq. (|3.16|) are determined from the 
initial conditions ( |3.15| ). 



In the general case the roots fcj are quite complicated, but in some special cases 
they reduce to a rather simple expression. For instance, in the absence of collisions 
one has IV = ^wi = 0> so that ec l- C P- 13| ) can easily be integrated. The probability 
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eq. ( p,12|) takes in this case the form: 

P(t) = -exp(-n)[cosh( V / T^A^t) + j^=^ sinh(y / r 2 - A%j]} 



(3.19) 

This result, which here followed from the general kinetic approach, was previously 
derived in [Tl| by using a different method. One can also readily see from eq. ( |3.19 ) 



that in the absence of a magnetic field our approach gives rise to the standard MSW 
result for the neutrino conversion probability P^fl : 

P(t) = ^ sin 2 ^ . (3.20) 

m 

4. Shrinkage of the flavour spin 
4.1. Solution for B = 

Let us first apply our general formalism to the simple case of vanishing magnetic 
field by setting T = 0. This has been previously studied in the literature by using 
the one-particle approach. The new effect revealed by the quantum kinetic approach 
is the shrinkage of the flavour polarization vector, which results in a more stringent 
bound on the vacuum oscillation parameters than the previous ones. The effect was 



first discussed qualitatively by Stodolsky in p5| . 

We take A m ~ y( vec ) ^> Y w , i.e. we assume that the system is far from the 
MSW resonance. We then find 



Sl = 7^ [l + ^ + 2A^ GAT" + °^ w/ ^ 



S>2 



A m , r w V3 , r 2 , - r 2 , v w v 2 V27 



i w yo i Wl - i w i w v x yzi / 4 



We thus have Sl + s 2 = 2T w /3 - T W (V 2 /A 2 J and 8l - s 2 = (2/ v / 3)A m (l + (T 2 W1 
r 2 y )/2A 2 n ), where in the difference we can neglect the dispersion term (D 2 ) — (D) 2 
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Ty V1 — Tyy- Substituting s± and s 2 into the roots ki given by eq. ( p.!8| ), we obtain 
ki = — T W (V X / A m ) 2 , k 2 = —T w + iA m and k 3 = —Y w — iA m . The transition 
probability ( |3.16| ) is then given by 

1 V 2 
Pu e ^u s = -{l~exp(-(VjA m ) 2 T w t) + ^[exp(-(VjA m ) 2 T w t)-e- rwt cosA m t\} . 

m 

(4.1) 

For times much larger than the collision time T^/, but much less than the depolar- 
ization time 

U = , (4-2) 

neutrino harmonic oscillations die away and the probability for v e — > v s conversion 
is aperiodic: 



P=- 

2 



t 



t 



1 -exp ~ — . (4.3) 



2t d 



This probability is small, 



but nevertheless the aperiodic term dominates over the contributions from the last 
two terms in the brackets in eq. (|4.1|) . 

The earlier result for the transition probability [|TTJ reads, 



P = V 2 12 A 2 = Sin2 200 (4 5) 

VJ^rn 2(1-2XCOS2# + X 2 ) ' ^ 

where x is defined through 

MeV = ( 2eV 2 ) " ^ 

As can be seen, our result eq. ( |4.3| ) differs from eq. (|4.5| ) by an additional large 
factor T w t ~T W /H ~ (T/MeV) 3 . This is due to spin shrinkage which is manifest 
in the kinetic approach. 
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4.2. Heuristic derivation of spin shrinkage 

The result eq. ( |4.3|) has been discussed earlier by Stodolsky fl25| , who also pointed 



out that a large conversion probability P = 1/2 may never be reached because de- 
polarization time td could be very large, in our case larger than the age of universe. 
Depolarization time can be estimated in a simple manner by noting that flavour spin 
rotation turns the longitudinal part of the spin into transversal, and at each collision 
the transversal part is, in effect, wiped out. This results in a shrinkage of the spin 
vector, and as t — > oo, P Ve ^ Va (t) approaches the value 1/2. The length of the spin 
vector P gives the degree of coherence, and as t — > oo, P vanishes, which corresponds 
to a completely mixed incoherent state. 

Qualitatively one can see this by writing the QKE (|3.3| )- (|3.6| ) in the familiar 
form p5| 



dP 

— = V x P - DP ± , 

dt 

where now V = V^ vec ^h z + V x h x , D ~ T\y and the transversal spin is given by 
P_l ~ P{V X /V Z ) with V x = Asin20 o . It then follows that 

dP 2 2 
— - -2T W P ± , 

yielding a shrinkage rate P~ 2 dP 2 /dt = —Yw x V 2 / (V^ vec ^) 2 , in agreement with eq. 
( [4.31) . Thus, in the presence of collisions the one-particle matrix density 

Pu(t,k) = ^[P (t) + P i (t)a i ]f v (k) 

tends in the formal limit t — > oo to the unpolarized one because of the slow spiralling 
of the spin vector into the origin: 

Pu(t,k)^^P (oo)f v (k) . 

The outcome is a unique equilibrium point for the system, given by P x = P y = P z = 0, 
where the number of active neutrinos equals the number of sterile neutrinos. 
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4.3. Big bang nucleosynthesis constraint 

The aperiodic term modifies considerably the constraint on the neutrino mixing pa- 
rameters. Let us look at this matter more closely by first recapitulating the derivation 



of the earlier result [ TT , and then presenting a similar deriva tion within our present 
approach. 

The condition for non-equilibration of sterile neutrinos is T s = PTw < H, which 
implies the constraint 

(^) = 0.6 sin 2 29, x y/W \ Am? \ /(2 eV 2 )/^) < 1 , (4-7) 

\ xi / max v 

where y = ±sc depending on the sign of the mass difference Am 2 = m 2 — m 2 .. The 



function f%(y), defined as 



= 1 ^9 ^00 I 2 > ( 4 - 8 ) 

1 + 2y cos 20 Q + y 2 



is the product of the probability (|4.5|) and the ratio Tw/H. For a small mixing angle 
it reaches the maximum value f^ ax ~ 3-^3/16 at y = (y/3 + cos 2 29 Q - cos2fl )/3. 



Substituting these values into eq. (|4.7|) one obtains the following allowed region for 



the mixing parameters [11]: 



| Am 2 | sin 4 2^o < 5 x 10~ 6 eV 2 . (4.9) 

For a typical mixing angle sin 29 ~ 0.1 this yields | Am 2 | <5x 10~ 2 eV 2 . 

In our case, where the conversion probability is given by eq. (|4.4j) , one obtains 
the constraint 

(^) = 2.5 x 10 6 sin 2 29, x L^lf™* (y) < 1 , (4.10) 

V XI / max eV 



where the function 



^ = i.o V oft I 2 
1 + 2y cos 29 + y 2 



14 



has the maximum value f™ ax = 1/4 cos 2 8q at y = 1. Thus, in our case the allowed 
region is 

| Am 2 | sin 2 29 < 1.6 x 1(T 6 eV 2 . (4.12) 

One should notice the different sin 28q dependence compared with the earlier result 
( f4.9| ). This difference gives rise to a more stringent constraint on the neutrino mass 
in our case. For example, for sin 29q ~ 0.1 we find | Am 2 | < 1.6 x 10 -4 eV 2 . The 
old and the new constraints in the (Am 2 , sin 26 , )-plane are depicted in the Figure. 

5. Constraints in the presence of a magnetic field 
5.1. Conditions for large magnetic effects 

Let us now consider neutrino propagation in hot magnetized plasma as described by 
the master equation (|3.13| ). In this section we will assume a non- vanishing magnetic 
damping T in eq. fl3.13| ), together with a non- vanishing quantum damping Yw- We 
assume that the neutrino system has not yet reached the MSW resonance region, so 
that the oscillation frequency in matter A m can be approximated by A m ~ \/( vec ). 
Furthermore, we take A^ 3> T 2 ^ , T 2 V1 . We also assume that the dominant plasma 
effect on neutrino prapagation is due to random magnetic fields as given by eq. (|2.1| ) . 
This is the case where the magnetic damping rate is much larger than the oscilla- 
tion rate of neutrinos, T ^> A m , and simultaneously A^/TTw -C 1. With these 
approximations, we are able to solve the QKE ( |3.13|) analytically. 

Let us study more closely these approximations. We may note that the magnetic 
damping rate 

T = 3fi 2 eff B rms {L = I h )L /Att{3 - 2p) 
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3.5(2.3)* 12 ^ p BLs(L = Ih) ( T , 5+2p /So ] 

l>l wm - xlU I j MeV (5.1) 



has a minimum value 

7r(3-2p) 10 24 G VMeV^ V T o 

which corresponds to the minimum value of the magnetic domain size L™ n , given 
in fl2j). 

Let us now make a further assumption. To illustrate the potential importance 
of the magnetic effects, we take B ~ Tq, which we believe is the largest possible 
random field that can be supported by the plasma. An example of this kind of a 
situation is possibly provided by the electroweak phase transition, where it has been 
argued |T7| that fluctuating Higgs field gradients could generate magnetic fields of 
the order of T^ w . 

To see when magnetic effects are important, we compare magnetic damping with 
IV ^ 4.0G|T 5 and the oscillation frequency A m ~ V^l Since > T w , a 

sufficient condition for the domination of magnetic damping is T min > V^ vec ). It 
turns out that even more restrictive condition is the requirement A 2 m /2T m i n Tw <C 1 . 
This condition implies the following range of validity of the approximations discussed 
above: 



i 
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T > (Mi^Mil x 10 -* . (5.2) 



MeV 2.3 

We note in passing that if we discard the hypothesis of the relic origin of the 
galactic seed field, thereby allowing considerably smaller domain sizes than given in 
eq. ( |2.2| ), the conditions T 3> A m ^> Tw, A 2 n /2TTw 1 can be fulfilled only for very 
small values of the index p, i.e. when the magnetic field is very close to a constant. 
Our analysis does not cover this case. 
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5.2. Application to primordial nucleosynthesis 

Assuming the validity of the approximations discussed in the previous subsection, we 
obtain for the roots (|3.18|) the following simple expressions: 

fci ^ -2T -Y w + A 2 m /2Y ~ -2T , 

k 2 ~ -T w - H ~ -r w , (5.3) 



I/ 2 A 4 V 2 

3 



fc = x + -m _ ' - 



2r i6r 2 iv 2r ' 

In k 3 we have neglected the second term A^/16r 2 rV as it is resonable to assume that 
V x < A m . From eq. (|5.3| ) one can then easily see that for t ^> only the last term 
in eq. (|3.12 ), the one proportional to C^e^, survives, so that the active-to-sterile 



neutrino conversion probability now takes the form 

P~i[l-e"^/ 2r ] . (5.4) 

In order to avoid a conflict with the primordial nucleosynthesis constraint, the pro- 
duction rate T s = PTw of the sterile neutrinos should be smaller than the expansion 
rate H ~ t^ 1 of the universe in the temperature range 1 MeV < T < Tqcd- When 
T < 1 MeV, left chiral neutrinos are decoupled, and collisions can no longer excite 
sterile neutrinos. Above the QCD phase transition temperature Tqcd sterile neutrino 
production is harmless, because their number densities are subsequently suppressed 
by the relative heating of the interacting species. From eq. ( |5.2| ) we find that, within 
our approximative solution, nucleosynthesis constraint is relevant only if p < 0.6. If 



this is so, we obtain from eq. ( |5.4j ) the condition | V x \ < 2HyT rnin /Tw- This can 
be translated into a constraint on the neutrino mixing parameters: 

1 Am2| sin2g <8xl(T 2 -^x ^ (-f-Y +P (p < 0.6) . (5.5) 



eV 2 u ~ v /^T^yV M eV> 

Since this limit is more restrictive the lower the temperature is, we substitute into 
eq. ( |5.5|) the lower bound of the validity range defined in eq. (|5.2| ). In this way we 
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obtain our final result 

I /\tH , I 1 33 21 

1 - ' sin 2^o < 2(3?r(3-2p))^ x 10"^. (5.6) 

This result is displayed for p = 0.5, 0.4 and 0.3, together with the B = case, in 
the Figure. 

We may note that the presence of a magnetic field implies a limit which is more 
stringent than if B = 0. Naturally, the limit is also more model dependent. In 
particular, it very much depends on our assumption that the primordial field is the 
seed field for the galactic dynamo, which reflects itself in T min . Moreover, our results 
are valid only for rather weak scaling behaviour. For larger p, it is likely that the 
bound ( |5.6|) would be much weaker. On the other hand, weak scaling is needed if 
primordial fields are to provide the seed field. 



6. Discussion and conclusions 

We have considered active-sterile neutrino conversions in hot plasma of the early 
universe. Our starting point is a general quantum kinetic equation, and in the col- 
lisionless limit we reproduce the results of The great advantage of the kinetic 



approach is that one is able to follow the evolution of the transition probability over 
any number of collisions. In particular, it allows us follow the chain of collisions over 
time scales of the order of the Hubble scale. This fact provided us with a new an- 
alytic solution which actually dominates the active-sterile conversion probability in 
isotropic plasma. It is related to the shrinkage of the flavour spin. Such behaviour 
is the combined effect of neutrino oscillations, which convert longitudinal spin com- 
ponents to transversal ones, and simultaneous damping of the transversal component 
by collisions. Qualitatively this effect, progressive decoherence, was discussed earlier 
in 1 25]. Here it was realized explicitly in the analytic solution ( |4.1|) . 
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An analogous shrinkage regime exists for Dirac neutrino spin-flip in plasma with 



a large scale random magnetic field [26] . The physical reason of the shrinkage of the 
total real spin there is the same as in our case. 

It is perhaps surprising that the inclusion of a magnetic field, which leads to an 
increase of the energy gap between active and sterile neutrino spectra (see eq. (|2.5|) ), 
tightens the constraint on the neutrino mixing parameters. In fact, the transition 
probability does decrease in the presence of a magnetic field exactly because of the 
energy gap increase. However, in the random field case the averaged conversion rates 
depend also on the neutrino squared mass difference Am 2 in a subtle way. While in 
the isotropic case there is a saturation of the conversion probability as a function of 
Am 2 in the presence of a magnetic field there is a linear dependence upon Am 2 . This 
behaviour is reflected in the stronger mixing angle dependence. 
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Figure Caption 

The nucleosynthesis constraint on the v e — v s mixing parameters, (a) old result; 
(b) new result, which includes the shrinkage of the spin vector. Shown are also the 
limits obtained in the presence of a primordial random magnetic field with the scaling 
parameter p = 0.3,0.4 and 0.5 (solid lines). In all cases the allowed region is below 
the contour. 
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